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Abstract. The Landau-Zener formula provides the probability of non-adiabatic 
transitions occuring when two energy levels are swept through an avoided crossing. 
The formula is derived here in a simple calculation that emphasizes the physics 
responsible for non-adiabatic population transfer. 



PACS numbers: 03.65.Ca, 31.50.Gh 



Energy level crossings have been of continued interest in physics, beginning 
with early analyses of the Born-Oppenheimer approximation 1 upto present day 
applications to the production of ultracold molecules [2] and the understanding of 
quantum phase transitions [3.. The probability of non-adiabatic population transfer 
during passage through an avoided level crossing is given by the well-known Landau- 
Zener formula (LZF) 4 . The usual derivations of the LZF involve the use of 
mathematical techniques such as contour integrals [4j [5] or Weber functions [1] , which 
might present a barrier to introducing non-adiabatic transitions in a basic quantum 
mechanics course. However the LZF can be understood using simple undergraduate- 
level physics, without getting mired in complicated mathematics. In the hope that 
it may be pedagogically useful and provide some physical insight into non-adiabatic 
population transfer, here is an elementary derivation of the LZF. 

Let |0) and |1) be two basis states of a system. The state vector is \tp) = 
V>o|0) + Assume that, at the initial time t = tj, the system is prepared in 

the |0) state: ">po(U) = l>^i(*i) = 0. Let the excited state |1) be higher in energy 
compared to |0) by an amount wo (we set ft = 1 everywhere for simplicity). Let these 
two states be coupled by a general time-dependent interaction. In the rotating wave 
approximation, the Hamiltonian for this two state system can be written as 

Here f2 is the strength of the time-dependent interaction, and u its frequency. 

Following Zener pQ , consider a process where the detuning A = to — loq is varied 
linearly with time: 4rA{t) = A = constant. This can be accomplished by varying the 
frequency u) of the perturbation fi e - *"*, or by tuning the energy level splitting u) 
(using a magnetic field, for example). For a constant detuning A, the dynamics of the 
system can be completely described by the usual Rabi flopping formula for a 2-state 
system [B] , but the Rabi formula is not a solution of the Schrodinger equation when the 
parameters ljq, lu are themselves functions of time. The exact analytic solution when 
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Figure 1. Energies of the dressed |0},|1) states in the rotating wave 
approximation. When the detuning is swept linearly through resonance, the 
probability of jumps across the avoided crossing is given by the Landau-Zener 
formula. 



A is linearly ramped in time involves Weber functions [1 . However, we shall find that 
the physics of this process can be understood without recourse to the mathematics of 
special functions. 

The main simplification arises because a time-varying detuning leads to rapid 
dephasing of coherent Rabi oscillations. Say we model the linear ramp as follows: to 
varies discretely in very small steps, each of which lasts for some small duration 5t. 
Within each of these time segments of duration St the detuning is a constant and the 
coherent evolution of the system is described by the usual Rabi flopping amplitudes. 
To obtain the state vector of the system for longer durations, the amplitudes from 
consecutive time-segments have to be added together. However the different detunings, 
and therefore different Rabi frequencies, in each segment mean that these amplitudes 
get added together with lots of different phases. Due to this dephasing introduced 
by the detuning ramp, the population dynamics of the system can be described quite 
well by summing the transition probabilities per unit time instead of the amplitudes. 
We estimate the dephasing time td as the duration to accumulate a phase difference 
of ~ 2ir rad between Rabi oscillations in consecutive segments: Td ~ y^- We will 
find that exact knowledge of this dephasing time is not necessary to get to the LZF. 

Under the action of the time-dependent Hamiltonian H, the probability per unit 
time for excitation out of the |0) state is the transition rate T, given by the familiar 
expression [6] 

T = n 2 — -2-5-rr . (2) 

A 2 +7 2 /4 v ' 

Here 7 is the decay rate of the Rabi oscillations, which we set equal to the inverse 
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dcphasing time: 7 = 1/td 
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Numerical calculation 
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Figure 2. Comparison between the analytical estimate in Equation Q and 
a numerical integration of the Hamiltonian in Equation (I}. The numerical 
calculation is for a ramp rate A = 3f2 2 , and does not include any spontaneous 
emission or other dissipative mechanisms. The oscillations seen in the numerical 
calculation are residual effects of coherent dynamics that are not captured by the 
Landau-Zener formula. 

Let Po(t) be the probability of being in the |0) state at a time t. After a time- 
segment of duration St, the probability Po(t + St) of remaining in the |0) state at the 
end depends on the transition rate T(t) within that segment. We have 

P (t + St) = [1 - St] Po(t) « e~ r « 5t P (t). (3) 

Therefore the probability Po{tf) of remaining in |0) over the entire duration of the 
process, from an initial time t\ to a final time tf, is 

rtt 



P (t f ) = exp 1 



T(t) dt 



exp 



exp 



exp • 



n 1 

20 2 



7 



A(i) 2 + 7 2 /4 
A(U) A 2 - 



dt 



(7/2) ; 



a t an 



)]}■ 



(4) 



7/2 / V 7 /2 

The mathematical limits A(fj) <C 7/2 A(t/) correspond to the physical case 
where the detuning is ramped from far below resonance, through the level crossing at 
A = 0, ending up far above resonance. In this limit, Equation (0} reduces to the LZF: 
Po(tf) = exp(— 27r57 2 /A). In principle, population transfer from |1) to |0) at a rate 
given by equation (2) must also be considered. It can be neglected in practice because: 
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(a) with the assumed initial conditions, no considerable population is accumulated in 
|1) until after the resonance is crossed (see figure [5]), and (b) from that point on, the 
detuning from resonance only continues to increase, making the population transfer 
back into |0) increasingly improbable. 

Note that the population in |0) is completely transferred over to |1) for an 
infinitesimally slow ramp through resonance (A — > => Po(tf) —> 0). This is an 
example of adiabatic passage [6]. For any nonzero ramp speed, the probability of a 
non-adiabatic transition is given by the Landau-Zener formula. 
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